Transcriptomics: Lecture 1
Biotech 7005/Bioinf 3000
Frontiers of Biotechnology: Bioinformatics and Systems Modelling
The University of Adelaide
Frontiers of Biotechnology: Bioinformatics and Systems Modelling
The University of Adelaide
Welcome To Country
I’d like to acknowledge the Kaurna people as the traditional owners and custodians of the land we know today as the Adelaide Plains, where I live & work.
I also acknowledge the deep feelings of attachment and relationship of the Kaurna people to their place.
I pay my respects to the cultural authority of Aboriginal and Torres Strait Islander peoples from other areas of Australia, and pay our respects to Elders past, present and emerging, and acknowledge any Aboriginal Australians who may be with us today
Introduction To Transcriptomics
Introduction
- Postdoctoral Bioinformatician, Black Ochre Data Labs, Adelaide
- Working in collaboration with members of the SA Aboriginal community
- Multi-omics project to identify and address the underlying causes of high T2D rates and complications
- Using genomics, epigenomics, transcriptomics and other layers
- My focus is on the transcriptomics layer
Why Transcriptomics?
- DNA can be described as being like a giant book of instructions
- Some regions are defined as genes
- Originally considered to be the basic unit of inheritance
- Now commonly used to describe a region of DNA transcribed into RNA
By Thomas Shafee - Own work, CC BY 4.0, Wikimedia Link
What Is Transcription
Definition
Transcription is the process of making an RNA copy of a gene sequence
Steps of Transcription
- RNA polymerase binds to the promoter along with \(\geq1\) transcription factors
- RNA polymerase creates a transcription bubble
- separates the two DNA strands, breaking hydrogen bonds between complementary DNA nucleotides.
- RNA polymerase adds RNA nucleotides
- complementary to the antisense DNA strand.
- RNA sugar-phosphate backbone forms
- Hydrogen bonds of the RNA–DNA complex break freeing the newly synthesized RNA strand.
Eukaryotic mRNA Processing
- Nuclear mRNA have 5’ cap added
- Protects single-stranded mRNA from degradation
- Regulates nuclear export
- Promotes translation
- mRNAs are polyadenylated at the 3’ end
- Also protects from degradation
- Aids in transcription termination, export and translation
- Introns are spliced out as required
Alternate Transcripts and Isoforms
Transcriptome Resources
- Reference Transcriptomes & Genomes are now commonly available
- Incorporate experimentally derived & predicted sequences + loci
- Gencode2 provide highest quality for mouse & human
- Release 48 (GRCh38): 78,686 genes + 385,669 transcripts
- Other organisms from Ensembl, RefSeq, UCSC etc
- Zebrafish, Rat, Chicken, Drosophila, Wheat, Yeast, E. Coli etc
- Sometimes we build novel transcriptomes from specific tissues
- e.g. sea snake venom gland, shiraz fruit
Early Transcriptomics
Northern Blotting
- Northern blot (Alwine, Kemp, and Stark 1977) extended DNA-based methods (i.e Southern blot)
- Earliest single-gene method
- Gel Electrophoresis then hybridisation with labelled probe
- Requires some knowledge of RNA sequence
- Images scanned \(\rightarrow\) Densitometric Analysis for crude quantitation
- Possible for different isoforms to be detected
- Sequence dependent
RT-qPCR
The CT values as actually estimated to a decimal value
- “Gold-standard” for measurement of transcription levels
- Single gene \(\implies\) not a high-throughput technique
- Targets a single transcript region with specific primers to produce cDNA
\(\rightarrow\) Polymerase Chain Reaction (PCR) - Each PCR cycle approximately doubles the target region
- cDNA produced is identified using fluorophores
- Fluorescence doubles with each cycle
- Once fluoresence passes a detection threshold, the cycle number is recorded
- Known as the Cycle Threshold (CT) value
Sanger Sequencing
SAGE & CAGE
- First high-throughput quantification method was Serial Analysis of Gene Expression (SAGE) (Velculescu et al. 1995)
- mRNA \(\rightarrow\) cDNA using biotinylated primers
- cDNA bound to beads (using biotin) & cleaved
- 11mer “tags” were ligated into long sequenced using linker sequences
- Sequenced using Sanger Sequencing
- Deconvolution & counting
Microarray Technology
- My search last week showed 69000 public miroarray datasets in the GEO database
- I reviewed a Scientific Reports submission using public array data last month
- Microarrays represent the birth of modern transcriptomics
- Thousands of genes could be measured simultaneously!!!
- Tens of thousands of public datasets \(\implies\) still being mined
- Established concurrently with the Human Genome Project (1990-2003)
- Databases & complete reference sequences become widely available
- All require fluorescently labelled cDNA copies of RNA
- Hybridised to the array using probes for known sequences
- \(\uparrow\) fluorescence \(\implies \uparrow\) RNA abundance
Two Colour Arrays
- Two colour microarrays were printed microscope slides
- Known probe sequences were printed to the surface in defined locations
- 60-75mer oligonucleotide probes
- Highly customisable by project
- Two samples per array
- One sample labelled with Cy5 (Red)
- Other sample labelled with Cy3 (Green)
- Scanned at 570nm (Cy3) and 670nm (Cy5)
Single Channel Arrays
- Affymetrix Arrays became dominant
- Factory manufactured
- Standardised layout for each organism
- Single sample per array
- Only scanned at one frequency
- More genes/array
- 25mer probes targeting 3’ end of transcript
- Captured only intact transcripts
3’ Arrays
- Each 3’ exon targeted by 11 unique 25mer probes \(\implies\) a probeset
- Possible to detect different transcripts only if 3’ exons differ
- Perfect Match (PM) probes \(\implies\) exactly matches target sequence
- Known to capture off-target signal \(\implies\) non-specific binding
- Arrays included mismatch probes (MM) with a change at the 13th position
- Sometimes returned more signal than PM probes 🤪
- Gene-level expression estimate obtained taking a robust average across probeset
- Model also included BG signal (but not always from MM probes)
Whole Transcript Arrays
- Whole Transcript Arrays released by Affymetrix in mid-2000s
- Marketed as Exon Arrays and Gene Arrays
- Probes along entire transcript BUT \(\leq\) 4 probes/exon
Microarray Analysis
Single Channel Data
- We will have multiple arrays from each condition
- Biological Replicates (hopefully \(\geq4\) per condition)
- Want to find changed expression in response to our biological hypothesis
- Will the some arrays have higher/lower overall signal?
- Pipetting errors, hybridisation variability etc
- Two Initial Problems to solve
- Adjust for overall differences in signal \(\implies\) Normalisation
- Removal of Background Signal (non-specific binding + optical noise)
\(\implies\)Background Correction
Normalisation
- The variation here is technical
\(\implies\) not due to biology - Higher variance reduces power of statistical testing
- Can we reduce this?
- Quantile normalisation
References
Abugessaisa, Imad, Jordan A Ramilowski, Marina Lizio, Jesicca Severin, Akira Hasegawa, Jayson Harshbarger, Atsushi Kondo, et al. 2020. “FANTOM Enters 20th Year: Expansion of Transcriptomic Atlases and Functional Annotation of Non-Coding RNAs.” Nucleic Acids Research 49 (D1): D892–98. https://doi.org/10.1093/nar/gkaa1054.
Adams, Mark D., Jenny M. Kelley, Jeannine D. Gocayne, Mark Dubnick, Mihael H. Polymeropoulos, Hong Xiao, Carl R. Merril, et al. 1991. “Complementary DNA Sequencing: Expressed Sequence Tags and Human Genome Project.” Science 252 (5013): 1651–56. http://www.jstor.org/stable/2876333.
Alwine, J. C., D. J. Kemp, and G. R. Stark. 1977. “Method for detection of specific RNAs in agarose gels by transfer to diazobenzyloxymethyl-paper and hybridization with DNA probes.” Proc. Natl. Acad. Sci. U.S.A. 74 (12): 5350–54.
Bolstad, B. M., R. A. Irizarry, M. Astrand, and T. P. Speed. 2003. “A comparison of normalization methods for high density oligonucleotide array data based on variance and bias.” Bioinformatics 19 (2): 185–93.
Fang, Rui, Walter N Moss, Michael Rutenberg-Schoenberg, and Matthew D Simon. 2015. “Probing Xist RNA Structure in Cells Using Targeted Structure-Seq.” PLoS Genet. 11 (12): e1005668.
Shafee, Thomas, and Rohan Lowe. 2017. “Eukaryotic and Prokaryotic Gene Structure.” WikiJournal of Medicine, January. https://doi.org/10.15347/WJM/2017.002.
Shalon, D, S J Smith, and P O Brown. 1996. “A DNA Microarray System for Analyzing Complex DNA Samples Using Two-Color Fluorescent Probe Hybridization.” Genome Research 6 (7): 639–45. https://doi.org/10.1101/gr.6.7.639.
Velculescu, V. E., L. Zhang, B. Vogelstein, and K. W. Kinzler. 1995. “Serial analysis of gene expression.” Science 270 (5235): 484–87.
Wang, Zhong, Mark Gerstein, and Michael Snyder. 2009. “RNA-Seq: A Revolutionary Tool for Transcriptomics.” Nat. Rev. Genet. 10 (1): 57–63.